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In the two Higgs doublet model, there is the possibility that the vacuum where the universe resides 
in is metastable. We present the tree-level bounds on the scalar potential parameters which have 
to be obeyed to prevent that situation. Analytical expressions for those bounds are shown for the 
most used potential, that with a softly broken Z2 symmetry. The impact of those bounds on the 
model's phenomenology is discussed in detail, as well as the importance of the current LHC results 
in determining whether the vacuum we hve in is or is not stable. We demonstrate how the vacuum 
stabihty bounds can be obtained for the most generic CP-conserving potential, and provide a simple 
method to implement them. 
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The two Higgs doublet model (2HDM) [T] is one of the simplest extensions of the Standard Model (SM) of particle 
physics. It has a rich phenomenology, allowing for the possibility of spontaneous CP breaking, as a possible explanation 
for the matter-antimatter asymmetry in the universe. It has a richer scalar content, with two CP-even scalar particles, 
a pseudoscalar and a pair of charged scalars; possible dark mater candidates; and many other interesting features. 
For a recent review, see In light of the recent discovery at the LHC of a particle that closely resembles the SM 
Higgs boson [311], but which seems to show some deviations from its expected behaviour, we can finally use the 
experimental data to choose from the plethora of proposed SM extensions. In particular, the 2HDM has shown to be 
quite capable of reproducing the available experimental results [SHE]. 

The price to pay for the rich phenomenology of the 2HDM is a scalar potential which is much more complex than the 
SM's, and which possesses a greater number of free parameters. Although unknown, those parameters are not wholly 
unconstrained. For instance, to ensure the existence of a minimum in the theory, the 2HDM scalar potential needs 
to be bounded from below. This severely constrains the quartic scalar couplings of the theory |9]. It is also usually 
required that all amplitudes involving scalars preserve perturbative unitarity [TOl [H] - an analogue, for the 2HDM, of 
the Quigg-Thacker constraints [T21[T3]. These bounds, once again, strongly restrict the potential's parameters. Known 
experimental evidence also comes into play: electroweak precision data is used, via the S, T and U observables [T4l - 
118] ■ to impose bounds on the 2HDM parameter space; and measurements from B-physics experiments impose serious 
restrictions on the scalar-fermion couplings of the 2HDM. 

In the present work we will address a new class of bounds on the 2HDM scalar potential, related to the existence 
of a metastable neutral vacuum. The vacuum structure of the 2HDM is much richer than the SM's. In fact, for 
certain choices of parameters, it is possible to have a vacuum which spontaneously breaks CP invariance - indeed, 
that is the reason Lee first proposed the model in 1973. For other regions of parameter space the vacuum may break 
the electromagnetic U{1) symmetry - such vacua are to be avoided at all costs. And of course, for large regions 
of parameter values, the vacuum of the 2HDM is "normal" and breaks electroweak gauge symmetry but preserves 
both electromagnetic and CP symmetries. The 2HDM has however a final surprise in store: it is possible for the 
scalar potential to display two such "normal" minima, both of them breaking exactly the same symmetries. Then the 
possibility arises that one of those minima is the one where we currently live, where the scalars' vacuum expectation 
values (vevs) give elementary particles their known masses; but in the second, deeper, minimum the vevs are such 
that all particle masses are completely different. "Our" vacuum is then but a metastable one, and the state of lowest 
energy lies below us. We call this situation a "panic" vacuum. In a recent work (TH] we studied under what conditions 
this occurred, for a specific version of the 2HDM, namely a softly broken Peccei-Quinn potential f20'. The conclusion 
reached therein was remarkable: the current LHC results enable us to conclude that, for this version of the 2HDM, 
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our vacuum is not metastable, i.e. it is the potential's global minimum. And this regardless of any cosmological 
considerations. 

In the present work we will present in detail the conditions under which the most general CP-conserving 2HDM 
potential can develop two minima. And also what bounds one can impose to prevent the occurrence of a panic 
vacuum. We will analyse in detail the possibility of panic vacua in the most popular version of the 2HDM, namely the 
one where a Z2 symmetry has been imposed on the lagrangian, but that symmetry is softly broken. We will present 
the exact analytical expressions for the bounds one needs to impose to ensure absolute stability of the vacuum for 
this model. We will show what the current LHC results have to say about the status of "our" minimum in the Z2 
potential. 



I. THE VACUUM STRUCTURE OF THE 2HDM 



The most general renormalizable 2HDM scalar potential is written as 
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where the coefficients 771^2, ^5,6,7 can be complex. Whereas in the SM there is only one possible type of vacuum - 
which preserves both CP and the electromagnetic U{1) gauge symmetry, but breaks SU{2)l x J7(l)y, and which we 
call the normal vacuum - in the 2HDM there are three. The normal vacuum corresponds to both doublets acquiring 
real and neutral vevs, 
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which can, without loss of generality, be taken to be both positive. In the charge breaking (CB) vacuum the dublets 
have vevs given by 
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where all the Ci are real. Finally, in vacua which spontaneously break CP, the fields' vevs have a relative complex 
phase. 
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A priori, all of these different vacua could coexist in the potential, raising the possibility of tunneling between different 
minima. However, in |21[ I22j . it was shown that this is impossible. If Vjv is the value of the potential at a normal 
stationary point, and Vcb its value at a charge breaking one, it is possible to show that the difference of the potential 
depths is given by 
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where — v\+V2 and rii?fj± is the square of the charged scalar mass, both of these quantities computed at the normal 
stationary point. And if that stationary point is in fact a minimum, then m^-t > and as such Vcb ~ Vn > - the 
normal minimum is guaranteed to be deeper than the CB stationary point. In |211 122j it was further proven that the 
existence of the normal minimum implies that the CB extremum is necessarily a saddle point. 

A similar result is valid for the comparison of a normal and CP stationary points: the difference in depths of the 
potentials is given by [211 ES] 
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with m\ being the pseudoscalar mass at the normal stationary point. Existence of a normal minimum thus automati- 
cally gives Vcp — Vn > - again, the normal minimum is the deepest. And again, in this situation the CP stationary 
point is a saddle point ^. 

These results can best be summarized as a simple theorem: no minima of different natures can coexist in the 2HDM. 
In other words, if the reader chose a set of 2HDM parameters such that a normal minimum exists, there is no need 
to worry about the existence of a deeper CB or CP minimum. 

The keen reader will notice that the theorem mentions minima of different natures, which raises the question of 
knowing how many minima of each type can exist in the 2HDM. For the CB and CP cases, the answer is simple: for 
a given set of 2HDM parameters, the minimization conditions admit only one CB vacuum of the form of eq. ([3|, and 
only one CP vacuum of the form of eq. Q ^. 

But in certain situations, the minimization conditions allow for several non-equivalent normal stationary points. 
And it was shown [231 124j that in fact two of those solutions can be minima. In other words, other than the normal 
vacuum with vevs given by eq. ([2|, for which one as v\+V2 — (246 GeV)^, there exists a second normal minimum N' , 
with different vevs {vi, W2}- For this second minimum, the sum of the squared vevs takes a different value, smaller or 
larger than (246 GeV)^. And the two minima are not degenerate, in fact they verify E5] 



Vn' - Vn 



N 



N' 



{V1V2 - V2V[) 



I \2 



(7) 



where the quantity (to^±/u^) is evaluated at both minima. A'' and N' . This raises the possibility that our minimum, 
with V = 246 GeV, is not the deepest one. And in fact, for certain regions of the 2IIDM potential. A''' is found to be 
the global minimum of the model - a minimum where the exact same symmetries have been broken, but where all 
elementary particles have different masses. In that situation our universe could tunnel to this deeper minimum, with 
obvious catastrophic consequences. We call this situation the panic vacuum. 

Before we proceed, let us clarify a common misunderstanding about these coexisting minima. They are sometimes 
dismissed out of hand as irrelevant for our understanding of the 2IIDM, since obviously model-makers can simply 
choose to start in the global minimum of the theory and develop perturbation theory from that point. Well, that is 
not correct: we do not have the freedom to decide in which of both minima we are currently at. The point is the 
following: assume we have access to the most precise experimental data, from some future collider; this data shows, 
without margin for doubt, that there are two CP-even scalars, a pseudoscalar and a charged one - in other words, the 
2IIDM describes particle physics. Further, let us assume that, from the data, we can obtain all necessary information 
to reconstruct, precisely, the potential (we will discuss this in further detail in section |ll]) . At this point we have 
the complete potential and can look at the minimization equations. It is only then, after the potential's parameters 
are locked, that we can verify whether or not there is a second minimum, and if so if it is deeper than ours, and we 
are in the panic vacuum. In other words, we cannot decide to choose the potential's parameters such that ours is the 
global minimum; it will be the experimental data which will provide us with that information. 

With that clarification out of the way, let us go back to the coexisting normal minima, and their difference in depths 
given by eq. ([t]). Should we worry about the existence of a deeper minimum than ours? When this issue arose in 
supersymmetry, concerning dangerous charge and colour breaking vacua [27], the existence of the deeper minimum 
was only considered problematic if the tunneling time from our minimum to it was found to be inferior to the age 
of the universe - only in that case should one exclude the parameters of the theory which originate both minima as 
dangerous. In this paper we will present, in section|v] an estimate of the tunneling times between minima (a daunting 
task at best, even for the 2IIDM). And as we will show, in many cases the current LHC results are, remarkably, enough 
to exclude the existence of panic vacua - thus curtailing the need to compute any tunneling times, anyway. 

Another natural question one might ask concerns the compatibility of the metastable vacuum with the thermal 
evolution of the Universe. Is it natural - or possible at all - that the early hot Universe could end up in a metastable 
vacuum after cooling down from electroweak temperatures Tew'^ Indeed, the thermal fluctuations omnipresent at 
those high temperatures would preclude formation of a long-lived region which was not in the true vacuum. This 
seemingly bars the Universe from getting stuck in a metastable state. However, in models with a sufficiently complex 
vacuum landscape - including the 2IIDM - this description is not fully accurate. Temperature corrections to the 
free-energy density of the scalar field can be such that the relative depth of two coexisting minima changes its sign 
at a certain critical temperature Tc significantly below Tew- In the particular case of the 2HDM, this possibility 



^ Of course, it is possible to choose the values of the parameters of the potential to obtain a CP minimum - but the results we are 

discussing here imply that for those parameter values no normal minima can ever be found. 
^ Two solutions of the minimization conditions which are related to one another by gauge transformations are degenerate and as such 

taken as a single vacuum. Likewise, any other exact symmetry of the potential gives rise to the same situation. 
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was mentioned and investigated in [55]. That means that, when cooHng down from Tew, the Universe goes through 
electroweak symmetry breaking and then stays in the global minimum until T drops below Tc- After that, the Universe 
is in a metastable state, but the temperature is already too low to activate the thermal transition to the true vacuum 
and the tunneling rate is also too weak. This mechanism could be the origin of the panic vacuum today. We can 
call it the "vacuum freeze-out" in analogy with the well-known freczc-out phenomenon for various particle species in 
cosmology. The only difference is that the origin of the anomalously slow dynamics which drives the system out of 
thermal equilibrium is not the Universe expansion, but rather the tiny dramatic tunneling rate. Checking which of 
the panic vacuum points we find below are compatible with the vacuum freeze-out is a separate issue, which is left for 
future work. For the purpose of the present paper, it is sufficient to stress that this mechanism is present in 2HDM 
and does not require any fine-tuning. 

Eq. ([7]) has one major drawback - it is written in terms of the vevs of both minima. That makes it quite cumbersome 
to deal with, if one wishes to know whether one's minimum is the global one of the potential. In fact, both {vi, V2} 
and {v[,v'2} are, in general, the solutions of two coupled cubic equations. The ideal situation would be to have a set 
of conditions that the potential should obey to prevent the occurrence of panic vacua, written in terms of quantities 
pertaining to the minimum N alone. This is in fact possible, based on the work of ref. "2'Sl 1241 [25] . There, the generic 
conditions that specificy the existence of two normal minima in the 2HDM were obtained, as well as the means to 
answer the question of whether a given minimum is the global one. Those methods were recently applied to a simple 
version of the 2HDM, the softly broken Peccei-Quinn potential. In ref. [19' we provided very simple conditions, trivial 
to implement, that ensure the non-existence of panic vacua. We will now show how those conditions can be obtained 
for more general models, and what conclusions about the potential's minimum one can extract from the current LHC 
data. 



The most used version of the 2HDM, in a vast array of theoretical and phenomenological applications, has a Z2 
symmetry, <I>i — >■ $1 and $2 ~^ —^2- This symmetry was initially introduced to prevent the occurrence of tree-level 
flavour-changing Higgs-mediated interactions [301 [31] j and it eliminates, in the potential of eq. Q, the parameters 
771^2; -^6 and A7. However, in order to allow for a vaster parameter space after unitarity constraints arc put in, the Z2 
symmetry is softly broken by the re-introduction of the quadratic term proportional to m^2- what follows we will 
consider the case where we have further required that the only source of CP violation in the model is that of the SM, 
i.e. an explicit violation of CP by the Yukawa terms, originating a complex CKM matrix. As such all coefficients of 
the potential are taken as real Requiring that the potential has a stationary point with vevs given by eq. ^ is 
tantamount to solving the following minimization equations. 



where we have defined A345 = A3 -I- A4 -I- A5 . In terms of the soft-breaking term mf 2 , the masses of the CP-even scalars 
- the lightest mh and heaviest mn -, the pseudoscalar mass tua, the charged Higgs mass mfi±, the mixing angle a of 
the CP-even mass matrix and the angle /3 defined as tan P — V2 /wi , we have the following expressions for the quartic 



^ There is however much interest in the complex 2HDM, in which we relax this assumption and consider complex values for and 
As [32133]. See also [34]. 



II. THE CASE OF THE SOFTLY BROKEN Z2 POTENTIAL 
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Now, all of the quantities which appear in the equations above are, in principle, possible to measure in experiments. 
The physical masses can be obtained by looking at invariant mass peaks. To establish which of the neutral scalars is 
A we need only look at their decays to ZZ or WW - A will not have those. The angles a and /3 can be obtained from 
combined measurements of the decays of h and other scalars to ZZ^ bb and r+r^. Finally, the soft breaking term 
TO^2 can be extracted, for instance, from a precision measurement oi h, H ^ 77. 

Of course, even if all of these scalars were discovered, the LHC almost certainly would not be able to provide 
enough precision for accurate determinations of all the Ai, but the point we wish to stress is this: collider experiments 
are a priori sufficient to determine all parameters in eqs. ^ and, from those, the quartic couplings A^. Using the 
minimisation conditions ([S]) we would then determine the values of nrf^ and 77122, ^^'^ parameters of the potential 
would be uniquely determined from experiments. 

At this point, in possession of all of the parameters of the potential, we can go back to the minimisation conditions ([s]) 
and try to solve them for different values of the vevs vi and W2- The soft breaking term 771^2 renders an analytical 
solution of these equations impossible. Of course, they allow for the trivial solution, both vevs equal to zero - a 
maximum of the potential. But it has been shown [531 HI] that they can lead to several solutions, and at most two 
non-degenerate minima. In fact, expressing the charged Higgs mass in terms of the parameters of the potential (see, 
for instance, eq. (204) of [2]), we can rewrite eq. ([?]) as 



where once again we see the crucial importance that the soft breaking term has - without it, the minima N and N' 
would be degenerate. 

And thus, there is the following tantalizing possibility: in the future, a precise determination of the parameters of 
the 2HDM potential leads us, by solving eqs. to determine that the model has more than one minimum; and, 
for some choices of parameters, that the minimum the universe currently resides in is not the global one - the panic 
vacuum we alluded to in the introduction. We are now going to provide the reader with a set of simple criteria to 
answer the following questions: under what conditions can there be two normal minima in the potential? Under what 
conditions is our vacuum, with v = 246 GeV, not the global one? We will now write those conditions, postponing 
their demonstration until section HVl 



The softly broken Z2 2HDM potential can have two normal minima if the two following conditions are met: 
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where we have defined 
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and also 

(14) 

These conditions are necessary and sufficient conditions for the existence of four ^ stationary points in the potential - 
but do not guarantee that two of those are minima (see appendix |b]). Nonetheless, only under these circunstances can 
the potential have a maximum of two normal minima [21] ■ Notice that these are trivial extensions of the conditions 
considered in [19j for the softly broken C/(l) model, which is a particular case of the Z2 case we are considering here 
(with A5 =0). Notice that they can be written only in terms of the potential's parameters, without any mention of 
a specific vacuum. 



In order to study the importance of the bounds of eqs. ( 11 ), ( 12 ), we have performed a vast scan over the parameter 
space of the 2HDM. We have taken my, = 125 GeV, 125 < tuh < 900 GeV, 90 < tua, mH± < 900 GeV, -7r/2 < a < 



7r/2, 1 < tan;3 < 40 and 
from eqs. ^) obey 
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< 900 GeV . We demanded that the quartic couplings of the potential (calculated 
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so that the scalar potential is bounded from below (in the "strong sense" as defined in ref. [3S]). We also required 
that these quartic couplings are such that they satisfy perturbative unitarity [101 E] and the electroweak precision 
constraints stemming from the S, T and U parameters [HHHj- Our simulation consists of 700000 "points", each one 
corresponding to a different combination of potential parameters. 




FIG. 1: Distribution, in the {x,y) plane, of the parameter space of the Z2 model. Blue (black) represents parameter choices 
for which the potential has a single minimum. In yellow (light grey) are the parameters for which the potential can have two 
neutral minima. The red (dark grey) curve is the astroid, inside which points with two minima have to lie. 



Eq. ( |12[ ) defines, in the (x,y) plane, a region of space delimited by an astroid. All the points with two minima will 
have to lie inside the astroid and have rriii + fc^ 77122 < 0, all other points will necessarily have just one minimum. In 
fig. [1] we show how a generic scan of the softly broken Z2 2HDM includes many points where two minima are possible, 
represented with the color yellow (light gray), as opposed to the points where the minimum is unique, painted blue 
(black). Though hard to see in fig.fll there are some blue points inside the astroid - those for which the condition ( 11 ) 



* In fact they are eight, but since the potential is invariant under <I>i — > —$1, ^2 —^2, four of them are degenerate with the other 
four. This is a manifestation of the U{1)y transformation, and it should not be taken into account. 



7 



is not satisfied. In total, the yellow region consists of over 140000 points. The existence of two neutral minima is not, 
therefore, a curiosity to be dismissed off hand - a full one fifth of the model's parameter space (after sensible cuts) 
does not have a single minimum. 



B. Existence of a panic vacuum 



We have defined the panic vacuum as the following situation: our vacuum, which is caracterized hy v = 246 GeV, 
is not the global minimum of the potential. The panic vacua are therefore a subset of the regions of parameter space 
for which there are two minima, and the conditions under which that occurs were written in the previous section. 

Remarkably, and again based on the work of refs. [SSlUllUn], it is possible to write extremely simple conditions 
for the existence of a panic vacuum. Let us define the "discriminant" 

D = ml^iml^ - fc2TO22)(tan P - k) (16) 

where tan/3 — vijvx as usual, and written, of course, in terms of the vevs of "our" vacuum. The existence of a panic 
vacuum is thus summarised in the following theorem: 

Our vacuum is the global minimum of the potential if and only if D > 0. (17) 

Therefore, if we only wish to make certain that we are in the global minimum of the potential, regardless of the 
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FIG. 2; Panic vacua in the {x, y) plane. Blue (black) represents points for which the potential has a single minimum. In yellow 
(light grey) points corresponding to two neutral minima. Superimposed on those are red (dark grey) points, which represent 
panic vacua. They are bound by the red (dark grey) astroid curve. 



number of those minima, requiring D > is a necessary and sufficient condition. It is not necessary to verify the 
conditions shown in the previous section. 

Again, the existence of a minimum of the potential deeper than ours is not just some curiosity or restricted to 
a very narrow region of parameter space: in the scan of 700.000 points of the softly broken Z2 model we discussed 
in the previous section, we found almost 126000 points which corresponded to the panic vacuum scenario! Their 
distribution in the {x,y) plane, with the variables defined in eqs. (13 1, is shown in fig- [2] the red (dark grey) points, 
now superimposed over those shown in fig. [Tj correspond to the panic vacua. Notice that only the left-hand side of 
the astroid contains panic vacua. This can be demonstrated using the expressions shown in section |IV| 



C. LHC results and the existence of panic vacua 



At the time of this writing, the LHC has no hints whatsoever of the existence of more than one scalar particle. 
Nonetheless, as we are about to show, the current data can already tell us a great deal about the nature of our 
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vacuum, and the existence, or lack thereof, of a deeper minimum in the 2HDM. 

First, some explanations concerning what we are comparing with experimental data: the Z2 symmetry which we 
impose on the 2HDM scalar potential has to be a symmetry of the whole lagrangian, otherwise we would be dealing 
with a non-renormalizable theory. Thus, that symmetry needs to be extended to the Yukawa sector, and there are a 
multitude of ways to do that. In this work we limit ourselves to the two most popular possibilities: Model I, in which 
only the doublet $2 couples to fermions; and Model II, in which $2 couples to all up-type quarks, and $1 couples 
to all other fermions. These two models have very different scalar-fermion couplings (see, for instance, table 2 of [Ij) 
and very different phenomenologies. 

The LHC data most relevant for Higgs physics at the moment are the ratios between observed rates of the Higgs 
boson decaying into certain particles and their expected SM values. Assuming that what is being observed is explained 
by the 2HDM, we define the said ratio for a given final decay state / of the lightest CP-even Higgs boson h to be 
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(18) 



including thus both production cross sections a and the branching ratios (BR) of the Higgs boson. Here, we are 
considering all possible Higgs production mechanisms, but current LHC results already allow us to distinguish, in 
certain cases, between some of those. Namely: 

• The gluon-gluon (gg) production mechanism, in which two gluons, one from each colliding proton, produce a 
Higgs boson via a triangle of quarks (mostly tops, with a small percentage of bottoms). Accordingly, we define 
a quantity, considering only the cross sections of the gg process: 

^2HDM ^pp ^gg^h) BR^"^^' {h 
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(19) 



• The vector boson fusion (VBF) mechanism, in which quarks inside the protons radiate electrowcak gauge bosons 
V = W, which "fuse" to become a Higgs bosons. The VBF rate is thus defined as 

.2HDM ^pp ^VV ^h) BR^"^^' {h 
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Finally, a word on experimental constraints: we already mentioned the electroweak precision data from LEP that 
we included in our scan of the 2HDM parameter space, via bounds of the S, T and U observables. Those only provide 
constraints on the scalar sector of the theory (assuming no extra generations of fermions are present). But there are 
plenty of data from _B-physics which provide constraints on the fermionic sector of the 2HDM, and which need to be 
taken into account. We have used the latest updated bounds from [36J ^ . These translate into a limitation on the 
tan/S-m^fi plane. 

In fig. [3] we show, for both models considered, the results we obtained for the rates of the light Higgs h into two 



Model I 



Model I 





However, we have not taken into account the Bs — fi~^ fJ.~ bound shown in that reference, since the analysis presented therein seems 
specific for the MSSM. Neither did we consider the B — > D^*''t~Ut anomaly observed by the BaBar collaboration, due to lack of an 
independent confirmation of it. 
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FIG. 3: Rzz as a function of R-/-/ for models I and 11. Green (light) points are the total of our simulation, red (dark) ones 
correspond to the existence of a panic vacuum. The solid (dashed) lines correspond to the Icr (2cr) experimental limits. 

Z bosons versus the rate of h into two photons. In green (light grey) we show all the points obtained in our scan, 
with the above constraints. In red (black) are the points for which a panic vacuum occurs. Please notice that the 
density of points is so great that there are plenty of green points scattered in the middle of the red ones. Fig. |3] is 
most relevant for showing that there are regions, in the plane R^^-Rzz, which are free from panic vacua. The solid 
and dashed lines shown in the plots correspond to conservative lu and 2a intervals on the combined values for Rzz 
and R^^, Rzz — 0-93 ± 0.28, R-^^ = 1.66 ± 0.33, which we took from ref. [55] , 

Notice that we are in the early days of Higgs experimental results, and as such many of the current numbers (such 
as the apparent excess in the two-photon rate) may well change considerably over time. The plots shown in fig. [3] are 
"invulnerable" to such likely changes, as future, more precise, restrictions on R^^ and Rzz can be imposed on them 
quite easily. Still, it is instructive to consider the current experimental bounds and see that, already, we can say much 
about the existence of a panic vacuum. As we see from fig. [3j for Model I the panic points are well away even from 
the 2(7 bands, which include some non-panic region as well. Not so for Model II, some of the panic region is included 
at 2(7. Of course, there are plenty of non-panic vacua choices of parameter space still allowed by the current data for 
Model II. As such, in these variables at least, both models seem capable of describing the current data, but that data 
does not exclude the possibility, in Model II, of our vacuum being metastable. 

What about other observables, for which there is already considerable information? The LHC collaborations have 
been able to measure - with considerable uncertainty - the ratio of two-photon Higgs events stemming from gluon- 
gluon production alone, and from the VBF mechanism alone. They are correlated, and we use the results from the 
ATLAS experiment, namely their la and 2a ellipses in the R^^l^-R}!^^^ plane [37]. The results of our simulations are 
shown in fig. [4] for both Model I and II. In these plots we see the situation is not as clear-cut as in the previous 

Model I Model II 




FIG. 4: R-^^ including only Higgs production via gluon-gluon fusion, and vector boson fusion, for models I and II. Colour code 
as in fig. |3] 

observables: in Model I we cannot exclude, at 2ct, the existence of panic vacua ^; and in Model II, even the la bands 
include panic vacua solutions. Notice that the ellipses contain plenty of green/light grey, non-panic points as well. 
And in these variables Model II agrees with the data at the la level, and as such describes the current data better 
than Model I. 

Finally, to conclude this brief comparison with experimental data, let us look at the tt rates, which have recently 
been measured by both LHC collaborations [38|. The current results are compatible with the expected SM value, 
ATLAS measured Rrr = 0.7 ± 0.7 and CMS, Rtt = 0.72 ± 0.52. Bearing these numbers in mind, as well as the ones 
presented above for the two-photon rates, we present what we have found for Rrr as a function of R^.j in fig. [5] The 
TT data (we represent the ATLAS bounds, since they are less restrictive), taken at face value, tells us that panic 
vacua are at least 2a disfavored in Model I, and the model agrees with the LHC results (at 2cr, barely, in R^^, at la 
in Rtt)- As for Model II, again at 1-a we notice plenty of panic vacua solutions not excluded by the data; but for 



Of course, the panic vacua points which now seem possible have been excluded by fig. 3 
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FIG. 5: i?TT as a function of _R^.y for models I and 11. Colour code as in fig. [3] 



much of Model IPs parameter space, we have agreement with the experimental resuhs at la, with or without panic 
vacua. 

In short: the current experimental data can already tell us a great deal about the stability of the vacuum in the 
2HDM. For instance, a measurement of Rzz and R^^ very close to 1, with sufficient precision, would exclude the 
possibility of panic vacua. Furthermore, this section shows that the same parameters which produce panic vacua 
do not correspond to some exotic, uninteresting corner of the model. They also predict values of observables which 
are not absurd and indeed can fall into the bounds of current experimental results. This, we would argue, is reason 
enough to take the existence of these panic vacua seriously in any phenomenological study of the 2HDM. 



III. PANIC VACUUM BOUNDS FOR THE GENERAL CP-CONSERVING POTENTIAL 



The conditions establishing the possible existence of two minima, and in particular of panic vacua, are amazingly 



simple and elegant, when written for the softly broken Z2 model, cf. eqs. (Ill, ( 12 ), ( 16 ) and the statement ( 17|. They 
are simpler still for the softly broken Peccei-Quinn potential, see ref. ^W[7 But for the most general CP-conserving 
2IIDM potential they cannot be written in a concise analytical manner, at least not in the usual notation. There is, 
however, a simple "recipe" which can very easily be implemented numerically when performing a study of the 2IIDM, 



and in this section we will give it in detail. We will present the demonstration of all bounds shown here in section IV 
We will leave out of this work any 2HDM potentials with explicit CP-breaking terms, since in those the discussion is 
even more difficult. Further, we will write the CP-conserving potential in a basis where all of its parameters are real. 
The existence of such a basis is guaranteed by explicit CP conservation liDJ. 

First, a brief discussion on notation: though writing the 2IIDM potential in terms of doublets, as in eq. ([I]), is 
extremely useful for many calculations {e.g., everything dealing with the fermion sector), in some instances a different 
notation - in which the potential is written in terms of gauge bilinear invariants - is crucial. For instance, the 
comparison of values of potentials at different vacua (eqs. ([s])- ([t])) is simple to obtain in the latter notation, but 
extremely hard in the former. Likewise, the conditions for existence of dual minima, or panic vacua, are far easier 
to establish in the bilinear formalism, which we now introduce. A remarkable feature of this notation is the fact 
that the 2HDM potential has a hidden Minkowski structure, when written in terms of gauge invariant bilinears. This 
formalism was developed in [53J [2J . Although similar Minkowskian notation was used in [33 HI] , those works did 
not fully exploit the freedom of non-unitary reparametrization transformations that the Minkowski formalism alludes 
to. The gauge invariant bilinears form a covariant 4- vector in a Minkowski space, (fi — 0, . . . 3), where we define 



I$2 + $^2*1) = -2 Re ($1$2) , 
ra = - - $^$2^ 



(21) 
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The allowed vectors r'^ = (ro,— r^) fill the forward lightcone LC+ defined by tq > 0, r^^r^ > 0. The apex of 
this cone corresponds to the electroweak symmetric vacuum, its surface corresponds to the neutral vacua, while its 
interior corresponds to charge-breaking vacua. Following the notation of [231 [23] and the conventions of eq. ([T]) , the 
CP-conserving scalar potential may be written as (with standard Minkowski space conventions) 



V = 



-A r^r" 
2 f^'^ ' 



where the 4-vector and the tensor A^^ are given by 



1 



= (Mo , M,) = --{ml^ + ml^), Re {ml^) , , -(m^^ _ m?,) 



(22) 



(23) 



(of course, Af^ = (ATq , -Mi)) and 



A„ 



/ Ai + A2 , , 
/ — ^ — + A3 A6 + A7 

A6 + A7 



V 





Ai — A2 



A4 + A5 


Afi — At 






A4 — A5 




Ai — A2 
2 

Ag — A7 


Ai + A2 



A, 



(24) 



As mentioned earlier, we are working in a basis where all parameters are real, which causes the appearance of several 
zeros in A and M . With the notation established, here are the preliminary steps required to verify whether or not 
the most general CP-conserving potential can have two neutral minima, and if one of them is a panic vacuum. 

• The first step in our "recipe" is the diagonalization of the tensor of the quartic couplings. Due to the Minkowski 
indices, this is achieved via a combination of rotations and Lorentz boosts . But there is a much simpler way, 
trivial to implement: define the matrix A = A^", which is obtained from eq. (24) by simply fiipping the sign of 



the three last columns. The 4x4 matrix A has eigenvalues Aa (a — 0, 1, 2, 3) determined by the usual equation, 

det(A - kal) ^ 0, (25) 

with eigenvectors V'^°'^ which satisfy (no sum in indices) 

Ay(") = A^y^''). (26) 

Solving for the eigenvectors and eigenvalues of A is trivially implemented within any numerical calculation 
package. Since the matrix A is not symmetric anymore, its eigenvalues and eigenvectors are in general complex. 

• The next step is to ensure that the potential is bounded from below. This means that the eigenvalues of A must 
obey the following conditions: 



All of them must be real. 
Ao > 0. 

Ao > {Ai , A2 , A3}. 



(27) 
(28) 
(29) 



The eigenvectors obtained in ([26j) are then real and can be normalized in such a way that one of them is time-like, 
the others space-like. Meaning, if the eigenvector corresponding to the largest eigenvalue Aq, obtained in (26) 
is given by V^'^'^ — (wooi ^lOi ^^20: ^^30)1 its overall normalization is such that, with our conventions, 



-"10 ^^20 

whereas, for the other three eigenvectors we must have 



30 



= 1 



(i)|2 _ 



''3i 



= -1. 



(30) 



(31) 



The diagonalization of A^j^ does not preserve the kinetic terms of the scalars. But it does not affect any details of calculations of vacua, 
or the value of the potential at vacua, so we need not worry. 
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• We now build a rotation matrix O, with the eigenvectors F^"-' serving as its columns. Which means, with the 
coefhcients v used in eqs. (30l and (31 1, Oat = Vab- This matrix O satisfies 



O-'AO ^ diag(Ao, Ai, A2, A3). 



(32) 



With this 4x4 matrix, we can obtain 4- vectors and r^, evaluated at the vacuum with neutral vevs vi/\/2 
and vil\Pl, in the basis where k^y is diagonal. In other words, the quantities M° and f" are obtained by 



Ml 


Vm3 j 



-Re fm? 
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12) 




(33) 



And thus Mq — Mi — — M*, etc. Since we began with the CP-conserving potential of eqs. (23 1 and (24), we 



are guaranteed to obtain M2 — f2 = Now in possession of the values of the eigenvalues Aq, A^; of the rotated 
quadratic coefficients Mq, M^; and of the rotated vevs fp, f^, the necessary conditions for existence of two neutral 
minima are very simple to write: 

If Mo > and ^/s^ + ^ff < 1, with 



Ml (Ao - A3) M3 (Ao - Ai) 

X = , y - 



Mo (A3 - Ai) 



Mo(A3 - Ai)' 



then the potential can have two neutral minima. 



(34) 



We emphasize that these are necessary conditions for existence of two neutral minima (see appendix B]) - although 
they are necessary and sufficient conditions for the existence of four normal stationary points. Remarkably, though, 
we have a necessary and sufficient condition to verify the global nature of our minimum - to know whether our {vi, V2} 
vacuum is the global minimum of the potential, we need only do the following: 

Let us define a discriminant D, given by 



D = MiM^hh] 



Our vacuum is the global minimum of the potential if and only if D > 0. (35) 

It is simple to verify that this procedure leads to the conditions laid out for the softly broken Z2 model in section [TT] 
Unfortunately, the diagonalization procedure explained above renders analytical expressions for the bounds unviable, 
in the case of the most general CP-conserving potential. But the "recipe" we provided in this section is quite easy to 
implement in a numerical study. 



IV. PANIC VACUUM BOUNDS: A DEMONSTRATION 



We will now demonstrate how the conditions for the panic vacua presented in the previous section are obtained. 
We are assuming scalar potentials which are stable in a strong sense (in the terminology of [3S]): that is, K^^r^^r'^ > 
everywhere on and in the forward lightcone (apart from the apex). Potentials stable in a weak sense cannot have 
two-minima configurations (29J, so we do not consider them. 

It was shown in [531 [51] that, for potentials stable in a strong sense, Ap,^ can be always diagonalized by an S'0(1, 3) 
transformation. This corresponds to the "recipe" we presented in the previous section. We therefore assume that we 
have performed that diagonalization, and the tensor A^;^ is written as 

A^^ = diag(Ao, -Ai, -A2, -A3) , (36) 



Remember, though, that CP breaking vacua, which would have r2 ^ 0, are excluded from the start, due to the existence of, at least, 
one normal minimum. 
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with the Aa coefficients satisfying the conditions of eqs. (27)- (29) so that the potential is bounded from below 



Likewise, the 4- vector of dimension-two coefficients is now given by = (Mq , Mi , , M3), calculated using eq. (33 ) 



Although this diagonalization does not preserve the form of the scalars' kinetic terms, for purposes of determining 
the extrema of the potential that is not problematic. With the potential written in terms of these parameters, the 
minimization problem is reduced to search for points lying on the surfaces of the forward lightcone, LC'^ , which 
minimize the potential. 

We remind the reader that we are interested only in neutral minima, and in particular in the possibility of two 
neutral CP conserving minima coexisting. Thus, we will not worry about the possibility of CP, or CB, minima, as they 
cannot exist when a normal minimum does. For instance, spontaneous CP violation could only occur if A2 > Ai, A3 
and the Ma obeyed the following condition: 



Ml 



< 



Mi 



(A2-Ai)2 (A2-A3)2 (A0-A2) 



(37) 



Therefore, our parameters have to be such that either A2 is not larger than Ai and A3 or the condition of eq. (37) 



is not verified. But in fact, it is simple to show that requiring the pseudoscalar squared mass to be positive(a 
necessary condition for a normal minimum to occur) implies that A2 is not the largest of the A;. We therefore assume 
that we are in that situation. 

Let us then go through the demonstration of the several conditions necessary for the existence of two normal minima 
- eqs. (34 1 - and the construction of the discriminants which establish how panic vacua occur - eq. (35). 



A. Condition Mo > 

Here we show that, if Mq < in the A^,y-diagonal basis, then the potential has only one non-zero stationary point, 
which must be the global minimum. The method we use is essentially the same as in [3S] and 



Finding neutral extrema of the potential (22) - that is, with values of r'^ restricted to the surface of lightcone, 
T^r/j^ = - benefits from using a Lagrange multiplier (. We introduce an auxiliary potential V = V — and 
minimize it, with respect to both the and ^. The minimization conditions thus become 

Af.ur" - M^ ^ Cr^ . (38) 
Using the explicit coefficients of the potential in the A^i^-diagonal frame , the minimization conditions become 

(Ao - 0^0 - Mo , (A, - On = . (39) 

Notice that, since the potential is CP conserving and no CP spontaneous breaking is being considered, M2 — and 
f2 = 0, always. Thus in these equations i = 1,3. This system has therefore three independent variables: C and two 
fi - the value of fg is then expressed as the positive square root of ff. 

Now, since fo is necessarily positive ^, the first equation in ( [39| implies that, if Mq < 0, the solution is found for a 
value of the Lagrange multiplier C, > Aq. Notice, also, that the condition t^q = ^ f? can be rewritten as 



i=l,3 





~Aoy 






-aJ 


\ Mo / 



This is a single algebraic equation of fourth order in C. Even without solving it, we can extract much information 
from it . 

Let us vary the value of C from Aq to infinity. We see that the variables (C — Aq) /{( — Ai) increase, in a monotonous 



manner, from zero to one. The expression on the left of (40) is therefore a monotonous function of ^, and it grows 
from zero to a maximum equal to (^ MD/Mq. 
We can therefore conclude that: 



If Mo < and Af^M^ > 0, then {Y, Mf)/Mi < 1, and the equation ([40]) has no solution in the region C > Ao. 
Thus, the potential has no non-trivial extremum. 

The only extremum - the global minimum - lies at = 0, and no electroweak breaking occurs. This situation 
is therefore physically uninteresting. 



Taking f > means we are excluding the trivial solution, all f/j, = 0. 
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• If Mo < and M^,M^^ < 0, then {J2 Mf)/M§ > 1, and the equation Q has only one solution at C > Aq. This 
is the only non-trivial stationary point of the potential, and it must be the global minimum, as it corresponds 
to y < 0. 

Thus no two minima can occur in this situation. 
Therefore, Mq > is a necessary condition for the 2HDM potential to have two minima. Since ro > per definition. 



eq. (39 1 implies that ^ < Aq. 



B. The astroid condition 

Here we show that, for the 2HDM scalar potential to have two normal minima, the values of the potential's 
parameters must be such that we are inside a region of space limited by the astroid curve defined in eq. (351). We 
will use the geometric approach to counting solutions of the minimization equations ( 40 ) developed in . Please 
remember that we are only interested in extrema with vevs without any relative phase, and thus r2 — 0. Since the 
existence of a normal minimum precludes a CP-breaking one, this means that a more generic analysis would discover 
a greater number of saddle points. But restricting ourselves to the r2 = case has no impact on the counting of 
possible normal minima. The analysis has a subtlety related to the ordering of the eigenvalues A^. Let us start with 
the case where Ai > A3 



Eq. (40) can be seen as defining an ellipse. In fact, if we define the variables mi — Mi/Mq and TO3 — AI^/Mq, the 



semiaxes of the ellipse will depend on and will be given by 

ai = |Ai-C|/(Ao-C) , - IA3 - CI/(Ao - C) ■ (41) 



In terms of these new variables eq. ( 40 ) thus becomes 



which is clearly the equation describing an ellipse in the (mi , m^) plane. Let us consider the family of ellipses which 
we obtain when we take all values of the Lagrange multiplier, —oo < < Aq, and count how many times this family 



of ellipses crosses a specific point (a, b) in the mi x plane. When an ellipse passes by that point, eq. (40) has 
a solution, which means that the potential has an extremum. Then, the number of times the ellipses pass over the 
point (a , b) will give us the number of non-trivial extrema of the potential. Also, it has been shown (for instance, 
in LSSJ that the larger ( is, the smaller the value of the potential. 
So, as C changes from —oo to Aq, we have: 

• In the limit C — > — oo, the ellipse is just the unit circle. This occurs because the semiaxes ai and 03 tend to one 
in this limit (solid blue line in fig.[6j left). 

• As C increases, the ellipse shrinks. It shrinks most quickly along the direction of the smallest semiaxis - and 



since Ai > A3, then, according to eq. (41), this means that the ellipse is contracting faster along the axis m3 
(dashed red line in fig. [6] left). 

• As C increases further, the ellipse on the (mi, m3)-plane shrinks even more. At C = A3 the ellipse collapses to a 
line segment (dot-dashed black line in fig. |6j left): 

m3 = 0, |mi|<mt-^i^. (43) 

Ao — A3 

• Notice that over the interval —00 < C < A3 the ellipses sweep exactly once all points inside the unit circle. For 
instance, there is only one ellipse which passes by the point (a, b) marked with an "X" in fig.|6j left. 

A single crossing means that in this region there is only one non-trivial solution of the minimization equations for 
this region of values of C. 

For A3 < C < Ai the situation is different: 



And also Ai > A2, since we are precluding the possibility of a CP-violating minimum. 
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FIG. 6: Example of the geometric evolution of the ellipse in eq. (40 1 with the value of the Lagrange multiplier. On the left, 
the range — cxd < < A3; on the right, Aa < < Ai. 



• For C > A3 the line segment again becomes an ellipse. When ( increase further, A3 < ^ < Ai, the ellipse shrinks 
along the mi axis and grows along axis (solid blue and dashed red lines in fig. [6] right). 

• At C = Ai the ellipses collapse to another line segment (dot-dashed black line in fig. [6j right): 



mi = , {m^l < TTig = 



Ai - A3 
Ao - Ai 



> to; 



(44) 



• During this evolution, the ellipses sweep a certain region in the (toi , m3)-plane, and each point inside this region 
is crossed exactly twice (see, for instance, the point marked with an "X" in fig.joj right). 

• It is possible to show (see Appendix |b]) that, of these two crossings, one is a saddle point and, if one imposes 
further conditions on the parameters, the other can be a minimum. 

Finally, when Ai < C < Ap, the ellipse grows infinitely (the semiaxes oi and 03 tend to infinity when ^ — >■ Aq), and 
the entire (toi , TO3)-plane is covered. But each point of the plane is crossed by an ellipse one single time. And since 
this last crossing corresponds to the largest value of C to yield an extremum, it corresponds to the global minimum 
of the potential. Please notice: regardless of whether the previous stationary points exist or not, this one is always 
present, and it is always the global minimum of the potential. 

We see, therefore, that we have the possibility of two minima - one may lie in the region where A3 < ^ < Ai (if 
we impose further constraints on the parameters of the potential), another (guaranteed) for values of C such that 
Ai < C < Ag. Had we considered the case where A3 > Ai, we would have similar conclusions: the figures obtained 
for that situation would be analogous to those in fig. [6) but rotated by 90°, and both minima would be found for 
values of the Lagrange multiplier in the ranges Ai < C < A3 and A3 < C < Aq. Notice that all we have been able 
to ensure, with the two conditions (11 1 and (35), is the existence of four normal stationary points, with one of them 



being certainly the global minimum of the potential. To make certain that one of the remaining stationary points is 
the second minimum we would need further conditions. 

Two neutral minima are therefore possible if, in the (toi , TO3)-plane, the values of the parameters of the potential 
are such that we are in a point inside the region of space covered by intersecting ellipses such as the two shown in 
fig. |6) right. The region covered by all of those ellipses is delimited by the astroid curve given in eq. ( |34[ ). To prove 
this, we apply the standard method of finding the envelope of a family of curves - the family of ellipses in question 
is generically described by an equation of the form F{mi , TO3 , C) = 0, and we also need to consider the tangent to 
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these ellipses at each point, determined by dF/dC, = 0. Explicitly, we need to solve 



F(mi , ms , C) 



1 



(Ai-C)^ 



1 dF 

3 Ik 



(A3 ^ (Ao - Q? 

mi 1 



(Ai-C)^ (A3-C)^ (Ao-C)^ 



= 



To solve these equations we define 



cos Lp 



Aq-C 

Ai-C 



siniy9 



Aq-C 

Aa-C 



so that eq. (45 1 is automatically satisfied, and eq. (46) becomes 

A3 - Ao 



Ao - Ai 2 



A, 



C 



sin = 0. 



From these two equations we determine the value of C, which, substituted in eqs. (47 1, yields 

Ao - A3 



COS^ (/3 



Al-A; 



.3 Ao-Ai 

mi , sm ip = — TO3 

A3 — Ai 



(45) 
(46) 

(47) 

(48) 

(49) 



Recalling that m.^ = Mi/ Mq and substituting these expressions in cos^ ip + sin^ </? = 1, we obtain the astroid curve 
which appears in our criterion for existence of two minima, eq. (Isil) Q.E.D. 



C. Conditions for panic vacua 

The necessary and sufficient conditions for the existence of panic vacua can be considered by themselves, circum- 
venting the need to verify whether or not the potential has two minima. 

The demonstration is extremely simple. Let us begin with the case Ai > {A2 , A3}. As we have seen in the previous 
section, in this situation the several possible stationary points obey the following relations: 

• The global minimum occurs for a value of the Lagrange multiplier, Cgi such that Cg > Ai > A3. 

• If another, local, minimum exists, it can only occur for a given value of the Lagrange multiplier, Cl, such that 
A3 < Cl < Ai. 

And this is all the information we require. Recalling the minimization conditions of the potential, written in terms 
of the Lagrange multiplier (eqs. (39)), we define the discriminants. 



Di 



-hAh , D = . 



(50) 



Given the minimization conditions, eqs. (39), we can write 



Di = -hMi = (C-Ai)f2 , 

D3 = -f3Af3 - (C-A3)f2 . (51) 
These discriminants can be computed for any minimum, i.e. for any given value of C- Then, we see that: 

• In the global minimum C = Cg- 

• Given that Cg > Ai > A3, we will have Di > and D3 > 0. 

• Thus, at the global minimum, D — D1D3 > 0. 

• If the second, local, minimum exists, it occurs for C = Cl- 



Notice that a sign difFerence in one of the denominators is irrelevant, due to the squares involved. 
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• Since A3 < < Ai, we will necessarily have £>i < and D3 > 0. 

• Thus, at the local minimum, D = D1D3 < 0. 

And so we see that the sign of Di discriminates between the local and the global minima, while the sign of D3 does 
not. If it happens that the potential has only one minimum, it will correspond to the case Di > 

Suppose we now have Ai > A3. The demonstration for this case is analogous to the one we have just given, with the 
following differences: the global minimum is now at C = Cg > A3 > Ai; the local minimum, if it exists, corresponds 
to a Lagrange multiplier Ai < < A3; at the global minimum we will have Di > and D3 > 0, at the local one 
Di > and D3 < 0. Thus, in this case, the sign of does discriminate between the local and the global minima, 
but the sign of Di does not. In any case, D = D1D3 is positive at the global minimum and negative at the local one. 

In conclusion, the product of Di and D3 is a quantity able to discriminate between the two normal minima: if we 
calculate it at a given minimum and find D = D1D3 > 0, that minimum is the global minimum of the potential; if 



D < 0, the minimum is local. Thus the conditions shown in (35) are proven 



V. LIFETIME OF THE METASTABLE VACUUM 

If particle physics is described by the 2HDM and D < 0, we are in a metastable minimum, and there is the 
possibility of tunneling to the true vacuum of the model. It can be argued, though, that the existence of these 
panic vacua is not sufhcient reason to exclude the parameters of the potential which produce them. In fact, if the 
tunneling time to the true vacuum is superior to the current age of the universe, the existence of a panic vacuum is 
completely acceptable. If, on the other hand, the tunneling time is inferior to the age of the universe, that region 
of parameter space ought to be excluded since it predicts a vacuum which would have already decayed, contrary to 
current experimental evidence. 

As such, it is interesting to try to estimate, for panic vacua, the decay width density, T/V, for a universe in such 
a metastable state. This is an application of the classic calculation by Coleman [42] for a potential with a single 
scalar field, very well discussed in the book by Rubakov |43) . Let us consider a potential with two minima, such that 
their difference in depths is given by e; there is a maximum between both minima, with height d relative to the local 
minimum. The value of the potential at the local minimum is taken to be zero. Then, one finds that 

l=Ae-^, (52) 

where A is a small prefactor, in comparison with the exponential of —B. The calculation of B, even for the simple 
case of a single scalar field, is very involved and requires a series of assumptions and approximations. Namely, it is 
assumed that: the height of the barrier S and the relative depth of the minima e satisfy 6/e ^ 0.06 [32]; that the 
only path linking both minima passes through a maximum of height S (which is a trivial assumption for a potential 
depending on a single field); and, while studying the expansion of the bubble corresponding to the true vacuum as 
the universe tunnels to it, one considers the so-called "thin wall" approximation, considering that the border between 
the regions of the universe lying in different vacua, as the bubble expands, has no thickness |44j . In this case. 



where A is the quartic coupling of the scalar potential. To have a rough estimate of what happens in the 2HDM 
case, we take A = max|Afc|, with k = 1, ... 5. But notice that, since in our potential we may have two minima and a 
multitude of saddle points and maxima, S is not necessarily the height of a maximum relative to the false minimum, 
but rather taken to be the height of the lowest stationary point - saddle point or maximum - relative to it. 



Considering the numeric factors in eq. (531, and the fact that A is limited in size by unitarity considerations, if the 

I becomes 
arge, and 



ratio (5/e is larger than about 1, the quantity B will be quite sizeable, and as such the decay width of eq. (|52 
extremely small - which means that the lifetime of the false vacuum, the inverse of F, becomes extremely 
no tunneling occurs during the lifetime of the universe. On the contrary, sets of parameters which predict S/e smaller 
than about 1 ought to be excluded from the model's parameter space, since they produce a tunneling time smaller 
than the age of the universe. 



Also, in the case where Mo < and there is a single minimum with (^q > Aq, discussed in section IV A both Di and D3 arc guaranteed 
to be positive, given that boundedness from below implies Aq > {Ai , A3}. 



18 



We have computed this estimate for the hfetime for ah our panic vacua. Using the equation (40), it is easy to 
discover ah possible normal stationary points, thus determining both e and S. We find that the vast majority of the 
panic vacua do not lead to vacuum lifetimes larger than the age of the universe. Quite the contrary, they correspond 
to very small values of S/e - either in potentials with minima of extremely different depths, or barriers of small height 
- which would mean that these vacua would have decayed long ago. If one believes this estimate, then, they must 
definitely be excluded. A very small percentage of points - less than 3% of the total of panic vacua found - seems, 
however, to have a false vacuum lifetime larger than the age of the universe, and as such ought to be retained. Our 
conclusions remain essentially unchanged, since the percentage of "safe panic vacua" is indeed very small. 

We must however point out several shortcomings in this lifetime estimate. In fact, the thin-wall approximation 
which leads to obtaining the estimate of eq. (53) breaks down if 6/e < 0.1 - and for most of our panic vacua 5/e 
is well below this value. As such, this estimate is problematic at best. Also, and perhaps even more serious, the 
estimates shown herein, which are widely used in the literature, are obtained in a model with a single scalar field. In 
the 2HDM, however, even if we use the gauge freedom to exclude the would-be Goldstone bosons, we are left with 
a potential with five real scalar fields. Now, to calculate the decay width density, we have to determine the bounce 
trajectory between both minima. This corresponds to a classical analysis of the scalar action in Euclidean space-time. 
In a scalar potential which is five-dimensional in its field content, the determination of the bounce is impossible to do 
analytically. Instead, we assumed that this trajectory corresponds, as we already explained, to the one that goes over 
the lowest intermediate saddle point, following the path of steepest descent. Although this seems like a reasonable 
approximation, it is known (43j that in multidimensional problems the bounce trajectory sometimes avoids the saddle 
point - in simpler terms, there is an easier path between both minima which avoids the intermediate saddle point. 
Thus, the fact that we found points which seem to be safe for tunneling is not too reassuring - there may well be 
another path between both minima which leads to much smaller lifetimes of the panic vacua. 

The remarkable thing, though, is that we can, in large measure, sidestep any cosmological considerations and simply 
look at what the LHC data tells us about the nature of the 2HDM vacuum. As we see for model I all panic vacua 
are excluded at 2cr, regardless of what their estimated lifetime might be. If however one chooses to believe in the 
traditional lifetime estimate for the false vacua, one will still find that most of the parameter space where panic vacua 
occur should indeed be excluded from phenomenological analysis, given that the lifetimes of those vacua would be far 
inferior to the age of the universe. 



VI. CONCLUSIONS 



The 2HDM scalar potential has a very rich vacuum structure, and the possibility of the coexistence of two normal 
minima in the tree-level is well-established. In this work we have derived the conditions that the parameters of 
the potential have to obey so that two normal minima may exist. We have also shown how one can build simple 
discriminants which allow us to conclude whether or not a given vacuum is the global minimum of the potential. 
These discriminants take on a particularly simple form for the most used 2HDM potential, the model with a softly 
broken Z2 symmetry, with explicit CP conservation. In this case, the vacuum of the model is global if and only if 

D = m\2{m\^~ k'^ml2){tanP - k) > 0, (54) 

with k = ^ A1/A2. 

We have performed a thorough scan of the parameters of the softly broken Zi model and shown that the occurrence 
of two normal minima is not confined to a non-interesting corner of parameter space: two normal minima occur 
very often in this model, if one does a blind choice of parameters. Further, we have seen that current LHC data 
disfavour at the 2tT level the possibility that, for Higgs Yukawa interactions of Model I, the minimum of the model 
which has v = 246 GeV and rrih — 125 GeV is not the global minimum of the potential. Thus, and even before any 
cosmological considerations, current particle physics experimental data already tells us a great deal about the nature 
of the 2HDM vacuum. However, for Model II, current LHC data cannot exclude the possibility that the model's 
vacuum is metastable, even at the 2a level. This would mean that the model has a deeper minimum, and the universe 
could, through tunneling, eventually reach that minimum. 

Following the standard calculations on the subject, we performed an estimation of the tunneling times involved in 
transition between vacua. To go beyond these simple estimates is extremely complex, and involves several assumptions 
which might be debatable. Taken at face value, however, our estimate shows that the vast majority of panic vacua we 
found would have lifetimes much smaller than the age of the universe. Such a situation can definitely be excluded on 
anthropic grounds, and the only panic vacuum points which we should really worry about are those with lifetimes of 
the order of the age of the Universe. But, as we discussed, our simple calculations might give an unreliable estimate 
of the vacuum lifetime, and improving them is a very complicated task. Thus, it is actually very difficult to establish 
which of the panic points can be truly worrying and which can already be disregarded due to cosmological reasons. 
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The strong feature of our result is that we can exchide these points right now, on the basis of the current LHC data, 
even without the need to improve the rehabihty of the hfetime calculations. Another, possibly interesting, line of 
research, considers the existence of two simultaneous minima - but in the situation where no panic vacua occur, and 
the second minimum lies above "ours" . What cosmological consequences might there be in such a situation? Could 
the universe have rested, for a short while, in the upper minimum, before cooling down to ours? What observable 
consequences would that have, if any? 

The calculations performed in this paper were, of course, all undertaken at tree level. The importance of radiative 
corrections to these results cannot, clearly, be dismissed. One may try to perform a Renormalization Group improve- 
ment of these results, by using the 2HDM /3-functions (see, for instance, |1S]) to run the values of the parameters of 
the potential with the renormalization scale up to the mass scale typical of the problem at hand, to curtail the need to 
compute the one-loop effective potential. However, please notice that both minima might have vevs of very different 
magnitude, and as such each minima would have a very different "typical scale". It isn't clear, then, up to what scale 
one should run the parameters, since comparing the value of the potential at both minima at different scales seems 
wrong. In fact, it has been argued, in the context of charge breaking minima in supersymmetric models |46 | 147 ) that 
the only correct procedure in the occurrence of two very distinct scales in two minima one wishes to compare is to 
compute the one-loop effective potential at both minima, coupled with an RG improvement of all parameters. Such 
a calculation is a gargantuan attempt. In this paper we limit ourselves to drawing attention to the problem that 
is already present at tree level - the existence of metastable, neutral vacua - and defer calculation of the impact of 
radiative corrections to future work. 

The panic vacua conditions presented in this work have a distinct practical advantage: they are extremely simple to 
apply in a numerical study of the 2HDM. In fact, they are almost as simple to apply as the usual bounded-from-below 
conditions, which are mandatory in any phenomenological study, and which involve the numerical determination of 
the parameters of the potential. In section III we presented a method which can be applied to more general (but still 
CP-conserving) versions of the 2IIDM potential to obtain the corresponding panic vacua conditions. 

But even if one were to obtain precise estimates of the tunneling times to panic vacua, and found them to be very 
superior to the age of the universe - which would remove any need to panic, really - the conditions we present in this 
work would still be interesting, in that they would allow us to determine, from particle physics experiments alone, a 
very interesting cosmological information: the true nature of the universe's vacuum. 
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Appendix A: Demonstration of panic vacuum bounds for the Z2 model 



The "recipe" we give in section |III| is completely general and it is a good exercise for the reader to apply if to the 
case of the softly-broken Z2 case and obtain the expressions shown in section [H] However, in that case, there is an 
alternative manner to reach the same results, and considerably simpler, and that is what we will now present. We 
start with the Higgs potential for this model, 



V = m2i$|$i+m22$^$2-m22 [$I$2+*^$l] +^Ai($l$i)2 + iA2($^$2)' 
+ A3($I$l)($^$2) + A4($I$2)(*J$l) + ^As [($1*2)' + (*|$l)' 



(Al) 



with all coefficients being real. The tensor is not diagonal in this basis: 



A„ 



/ ^i±^ + A3 




\ 2 




A4 + As 

A4 - As 




2 




A1+A2 



A3 ) 



(A2) 
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In order to diagonalize it, we need to equilibrate Ai and A2. This can be achieved by rescahng the doublets: 

1/8 



Upon this change, the potential becomes 



+ A3($I$l)($J$2) + A4($I$2)($^*l) + ^Ag [($1$2)' + 



In this basis, A^j^ is diagonal and has the following eigenvalues: 

Ao = ^(VAlA^ + Aa), Ai = -i(A4 + A5) , A2 = -^(A4-A5) 



A, 



^(-VA^ + As). 



(A3) 



(A4) 



(A5) 



Notice how the conditions of eqs. (27)- (291 lead to the well-known bounded-from-below constraints. In this basis, 
the components of vector Af^ take form: 



Mo 



1 



1 



(mfiq- + m^2'Z" ) , Mi = mi^ , Af2 = , M3 = -{mi^q 



2 2\ 

miiq ) 



(A6) 



If the neutral minimum point is parameterized as 



($?) = — cos/3, (<i>°) = -^sin/3, 



then the components of (r^) in the A^^^-diagonal basis are 



(A7) 



(fi) = -w^ sin/3cos/3 , {f'i) = —{q^s\T?fi-q ^ cos^ /3) , 



The two discriminants then become 



4 7 

IV „_2 



Di = V sin /3 cos /3 , I?3 = — — cos /3 mf^ 




tan^ P - 




(A8) 



(A9) 



Since we pay attention only to the signs of Z?i, one can also redefine them by removing factors which are guaranteed 
to be positive. With trivial manipulations one then arrives at eq. (16). 



Appendix B: Classifying stationary points 



Throughout this work we have used the method of Lagrange multipliers to determine the existence of several 
stationary points. However, to discover whether they are minima, maxima or saddle points, one must clearly look at 
the Hessian matrix computed at each of the stationary points. In terms of the Lagrange multiplier ^, the Hessian is 
easy to calculate, and is given by 



H = 



Ao-C 

C-Ai 

C-A2 

C-A3 



(Bl) 



Recalling the results shown in section IV B let us consider, for starters, the case where Ai > {A2 , A3}. We have 



said that the global minimum occurs for a value of the Lagrange multiplier Ai < Cg < Ag - and indeed, we see from 
the expression for the Hessian that with C — Cg 8.11 the elements in the diagonal of the matrix are positive, and as 
such this stationary point is guaranteed to be a minimum. Likewise, the first stationary point - that with a value of 
Q smaller than A3 - gives at least two negative values in the diagonal, and as such cannot be a minimum. 

It is rather trickier to verify if one of the other stationary points is a minimum, or a saddle point. In fact, it is not 
necessary that all of the diagonal elements of the Hessian matrix of (Bl) be positive for a stationary point to be a 



minimum. One has to require a milder condition, namely that H be positive definite in the subspace tangent to the 
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stationary point, in the space of the r^. A necessary and sufficient condition for the existence of a normal minimum 
becomes rather elaborate, and we will not present it here, since it is not necessary for the purposes of this paper. But 
we draw the attention of the reader to two very simple necessary conditions that must be obeyed so that one has a 
normal minimum, in the softly broken Z2 model. It is simple to show that the squared masses of the charged and 
pseudoscalar Higgses are given, in terms of the Lagrange multiplier C and in the basis where the A tensor is diagonal, 
by 

m^± = Cfc^«^ , TO^ = (C-Aa)^:^^^ . (B2) 

Thus, we see that any normal minimum must correspond to a value of a Lagrange multiplier which satisfies C > and 
C > A2. This latter condition, by the by, is the reason why we had to consider that A2 was not the largest of the A^. 
However, to specify that a given stationary point is a minimum, we would also need to look at the squared masses of 
the CP-even scalars and require that they be positive - and that condition cannot be cast into a simple form in terms 
of the parameters of the potential. 
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